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ABSTRACT 

We analyze MegaZ-LRG, a photometric-redshift catalogue of Luminous Red Galaxies 
(LRGs) based on the imaging data of the Sloan Digital Sky Survey (SDSS) 4th Data 
Release. MegaZ-LRG, presented in a companion paper, contains > 10 6 photometric 
redshifts derived with ANNz, an Artificial Neural Network method, constrained by a 
spectroscopic sub-sample of ~ 13,000 galaxies obtained by the 2dF-SDSS LRG and 
Quasar (2SLAQ) survey. The catalogue spans the redshift range 0.4 < z < 0.7 with an 
r.m.s. redshift error a z « 0.03(l+z), covering 5,914 deg 2 to map out a total cosmic vol- 
ume 2.5 h~ 3 Gpc 3 . In this study we use the most reliable 600,000 photometric redshifts 
to measure the large-scale structure using two methods: (1) a spherical harmonic anal- 
ysis in redshift slices, and (2) a direct re- const ruction of the spatial clustering pattern 
using Fourier techniques. We present the first cosmological parameter fits to galaxy 
angular power spectra from a photometric redshift survey. Combining the redshift 
slices with appropriate covariances, we determine best-fitting values for the matter 
density Q m and baryon density £!b of il m h — 0.195 ±0.023 and Slb/f^ m = 0.16 ±0.036 
(with the Hubble parameter h = 0.75 and scalar index of primordial fluctuations 
"■scalar = 1 held fixed). These results are in agreement with and independent of the 
latest studies of the Cosmic Microwave Background radiation, and their precision is 
comparable to analyses of contemporary spectroscopic-redshift surveys. We perform 
an extensive series of tests which conclude that our power spectrum measurements are 
robust against potential systematic photometric errors in the catalogue. We conclude 
that photometric-redshift surveys are competitive with spectroscopic surveys for mea- 
suring cosmological parameters in the simplest "vanilla" models. Future deep imaging 
surveys have great potential for further improvement, provided that systematic errors 
can be controlled. 
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1 INTRODUCTION 

The large-scale structure of the Universe is one of the most 
important probes of cosmology, encoding information about 
the fundamental parameters of the Universe and the physi- 
cal processes governing the development of cosmic structure 
and the formation of galaxies. In particular, the large-scale 
power spectrum of matter fluctuations depends on a differ- 
ent combination of cosmological parameters to the Cosmic 
Microwave Background anisotropy spectrum. The combina- 
tion of these two datasets provided the first compelling evi- 
dence for the current standard model of cosmology, in which 
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two-thirds of the present energy density of the Universe is 
resident in some form of "dark energy" which is driving the 
acceleration of the cosmic expansion rate (e.g. Efstathiou, 
Sutherland & Maddox 1990; Ostriker & Steinhardt 1995). 
This model was spectacularly confirmed by subsequent ob- 
servations of distant supernovae (e.g. Riess et al. 1998; Perl- 
mutter et al. 1999). Measurements of large-scale structure 
have continued to constitute an important component of cos- 
mological parameter fits (e.g. Efstathiou et al. 2002; Percival 
et al. 2002; Spergel et al. 2003; Pope et al. 2004; Tegmark 
et al. 2004b; Seljak et al. 2005; Sanchez et al. 2006). 

The pattern of cosmic structure has been delineated by 
a series of increasingly impressive galaxy redshift surveys. 
The current state-of-the-art is represented by the 2-degree 
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Field Galaxy Redshift Survey (2dFGRS; Colless et al. 2001) 
and the Sloan Digital Sky Survey (SDSS; York et al. 2000), 
which have provided exquisitely detailed maps of the lo- 
cal (z ~ 0.1) Universe (e.g. Percival et al. 2001; Cole et 
al. 2005; Tegmark et al. 2004a; Eisenstein et al. 2005). At 
higher redshifts (z ~ 1), surveys by the Deep Extragalac- 
tic Evolutionary Probe (DEEP2; Davis et al. 2003) and the 
VIRMOS-VLT Deep Survey (VVDS; Le Fevre et al. 2003) 
have produced impressive measurements of the clustering 
pattern over an area of a few square degrees (e.g. Coil et 
al. 2004; Le Fevre et al. 2005). But despite these efforts, 
today's most pressing cosmological problems - such as the 
investigation of the nature of dark energy, the discrimina- 
tion between competing models of inflation, and the mea- 
surement of neutrino masses - demand the construction of 
vastly larger spectroscopic surveys of millions of faint galax- 
ies spanning thousands of degrees at high redshift. 

The pattern of large-scale structure can also be mea- 
sured by imaging surveys via its projection on the sky. Al- 
though the lack of redshift information weakens the precision 
with which the clustering functions can be determined, this 
is partially compensated for by the ease with which large ar- 
eas of sky can be surveyed. For example, the APM Galaxy 
Survey, constructed by scanning photographic plates cover- 
ing several thousand square degrees of the southern galactic 
cap, was used in the early 1990s to produce the best then- 
existing measurements of the three-dimensional large-scale 
clustering pattern (Baugh & Efstathiou 1993). 

Moreover, the power of imaging surveys can be greatly 
enhanced by the use of photometric redshifts, which are es- 
timated from broadband galaxy colours rather than from 
spectra. These redshift estimates enable the catalogued 
galaxies to be divided into many quasi-independent radial 
slices, greatly improving the signal-to-noise ratio of the re- 
sulting clustering measurements. The utility of photometric 
redshifts is now well-established, with many successful tech- 
niques being employed. The r.m.s. precision a z with which 
galaxy redshifts may be determined varies with the method 
and filter set used, together with the galaxy type, magnitude 
and redshift, but at best is currently cro = o~ z /(l + z) ~ 0.03 
(e.g. COMBO-17; Wolf et al. 2003). Photometric redshifts 
have already been used to construct volume-limited samples 
of low-redshift galaxies, and measure their angular cluster- 
ing properties as a function of luminosity and rest-frame 
colour (e.g. Budavari et al. 2003). 

Certain sub-classes of galaxy provide especially ac- 
curate photometric redshifts, in particular Luminous Red 
Galaxies (LRGs), for which the optical colours change very 
rapidly with redshift owing to a significant spectral break 
at 4000A ( see Eisenstein et al. 2001; Padmanabhan et al. 
2005). Moreover, these galaxies inhabit massive dark matter 
haloes (Zehavi et al. 2005) and are consequently highly effi- 
cient (although biased) tracers of the underlying clustering 
pattern. In this sense, a sub-sample of LRGs should provide 
the optimal set of tracers of the large-scale structure. Al- 
though by constructing such a sub-sample we are excluding 
a large fraction of the overall galaxy population, this is im- 
material (as far as determining cosmological parameters is 
concerned) provided that the survey is sufficiently deep that 
the contribution of shot noise to the error in the clustering 
measurements is negligible. 

In this study we perform a clustering analysis on a large 



sample of ~ 10 6 Luminous Red Galaxies selected photo- 
metrically from the imaging component of the SDSS 4th 
Data Release (DR4), spanning the redshift range 0.4 < z < 
0.7. We believe that this investigation is timely for several 
reasons. Firstly, we aim to demonstrate the feasibility of 
extracting high-quality large-scale structure measurements 
from an imaging dataset employing photometric redshifts. 
Several forthcoming imaging surveys are planning to pro- 
duce accurate cosmological measurements in this manner - 
for example, the Kilo-Degree Survey (KIDS), the Dark En- 
ergy Survey (DES), and the Panoramic Survey Telescope 
(PanStarrs) - whereas next-generation spectroscopic sur- 
veys still require the development of challenging new in- 
strumentation such as the Wide-Field Multi-Object Spec- 
tograph (WFMOS) or the Square Kilometre Array (SKA). 
The cosmological parameter constraints resulting from fu- 
ture photometric redshift imaging surveys have been simu- 
lated by Seo & Eisenstein (2003); Amendola, Quercellini & 
Giallongo (2004); Dolney, Jain & Takada (2004); Blake & 
Bridle (2005) and Zhan et al. (2006). 

Secondly, the recent detection of baryon acoustic oscilla- 
tions in the galaxy clustering pattern (Eisenstein et al. 2005; 
Cole et al. 2005; Huetsi 2006a) has provided a powerful new 
cosmological probe. The existence of this preferred scale or 
standard ruler in the galaxy distribution will eventually per- 
mit accurate measurements of the properties of dark energy 
(Blake & Glazebrook 2003, Seo & Eisenstein 2003). Eisen- 
stein et al. (2005) and Huetsi (2006a) identified the preferred 
acoustic scale in the distribution of LRGs in the SDSS spec- 
troscopic sample, with mean redshift ~z ~ 0.35. However, the 
baryon oscillation scale can in principle be extracted from a 
photometric-redshift survey, and the completed SDSS imag- 
ing survey could produce a significant detection at an inter- 
estingly higher effective redshift of z ~ 0.55 (Blake & Bridle 
2005). 

Thirdly, this dataset provides a useful test case for the 
comparison of different methods for measuring the large- 
scale structure. In this study we consider two techniques: 
the determination of the (two-dimensional) angular power 
spectrum in redshift slices via a spherical harmonic analysis 
(Peebles 1973); and a direct re-construction of the (three- 
dimensional) spatial power spectrum using a Fourier anal- 
ysis which only retains large-scale radial components (Seo 
& Eisenstein 2003). The direct Fourier treatment may be 
preferable for the detection of features in the power spec- 
trum because it avoids unnecessary binning of data along the 
line-of-sight (which is equivalent to a convolution or smooth- 
ing of the measured power spectrum). However, the direct 
measurement of the spatial power spectrum requires an as- 
sumption of values for the cosmological parameters, thus 
parameter-fitting becomes a more complex iterative exer- 
cise. 

Finally, this survey provides an ideal demonstration of 
the latest techniques for determining photometric redshifts. 
We use ANNz, an artificial neural network approach (Firth, 
Lahav & Somerville 2003) which compares very well with 
other successful photometric redshift methods (see Collis- 
ter & Lahav 2004). The neural network technique requires a 
"training set" of spectroscopic redshift measurements for a 
sub-sample of the photometric catalogue. We use data from 
the 2dF-SDSS LRG and Quasar (2SLAQ) survey for this 
purpose (Cannon et al. 2006). The resulting photo- z cata- 
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logue of ~ 10 LRGs, which we have named MegaZ-LRG, is 
described by Collister et al. (2006). 

This paper is organized as follows. In Section 2 we de- 
scribe the construction of the LRG catalogue from publicly- 
available SDSS data, the definition of the angular selection 
function, and the determination and validation of the pho- 
tometric redshifts and error distribution. Section 3 details 
the measurement of the angular power spectrum in redshift 
slices using a spherical harmonic analysis, the fitting of cos- 
mological models and the resulting probability contours for 
the cosmological parameters. Section 4 contains an exten- 
sive series of tests for systematic photometric errors in the 
imaging data that may potentially influence these cluster- 
ing measurements. Section 5 presents a direct determination 
of the spatial power spectrum via a Fourier analysis and 
demonstrates consistency with the angular power spectrum 
results. We summarize our conclusions in Section 6. 

We note that a parallel investigation of a similar dataset 
was carried out by Padmanabhan et al. (2006) using a dif- 
ferent photometric-redshift method and power spectrum es- 
timation procedure. Where our results can be compared, the 
agreement is good. 



2 DATA 

We analyze galaxy clustering in the MegaZ-LRG database, 
a photometric-redshift catalogue of Luminous Red Galax- 
ies based on the imaging component of the SDSS 4th Data 
Release. The construction of this catalogue is described in 
detail by Collister et al. (2006), and we only provide a brief 
description here. 

2.1 Selection criteria 

We selected Luminous Red Galaxies at intermediate red- 
shifts (0.4 < z < 0.7) from the SDSS imaging database 
using a series of colour and magnitude cuts (Collister et al. 
2006). These cuts were designed to match the selection cri- 
teria of the 2dF-SDSS LRG and Quasar (2SLAQ) survey 
(Cannon et al. 2006), a spectroscopic follow-up to the SDSS 
imaging survey using the 2-degree Field (2dF) spectrograph 
at the Anglo-Australian Telescope. These 2SLAQ spectro- 
scopic redshifts were used to train and test the photometric 
redshift code, which we then applied to the entire set of 
LRGs selected from the SDSS imaging database. A similar 
procedure was carried out for earlier versions of the SDSS 
and 2SLAQ databases by Padmanabhan et al. (2005). 

The 2SLAQ survey obtained good-quality redshifts for 
about 13,000 objects in selected fields of the SDSS equa- 
torial stripe (at declination 8 £s 0°). Applying the same 
2SLAQ selection criteria to the entire SDSS DR4 imaging 
catalogue returned ~ 10 6 LRGs. The 2SLAQ survey demon- 
strated that these selection criteria are « 95% efficient in 
the identification of intermediate-redshift LRGs. The most 
significant contaminant, accounting for virtually all of the 
remaining ~ 5% of objects, is M-type stars. In the 2SLAQ 
survey these stars, which have very similar gri colours to the 
LRGs, are readily excluded using spectra. This is not pos- 
sible for a photometrically-selected sample, but star-galaxy 
separation cuts can be used to reduce the stellar fraction 
further, as described below. 



The 2SLAQ selection criteria fluctuated a little at the 
beginning of the survey. Specifically, the faint limit of the 
i-band magnitude idcv, and the minimum value of d pe rp (a 
colour variable used to select LRGs), were varied slightly. For 
the majority of the 2SLAQ survey, the criteria idev < 19.8 
and ciporp > 0.55 were used. In this paper we analyze a 
conservative version of the photometric-redshift catalogue 
in which these limits are strictly enforced, avoiding any ex- 
trapolation of the "training set" calibration. 

We note that our training sub-sample is extrapolated 
in sky position. The 2SLAQ targets lie exclusively in the 
equatorial stripe at declination 8 « 0°, so may not fully 
trace effects such as dust extinction which depend on sky 
co-ordinate. However, we find our results to be independent 
of varying dust extinction, as described in Section 4. 

We applied some star-galaxy separation criteria to the 
resulting catalogue (Collister et al. 2006). Firstly we made 
cuts based on galaxy angular size and the difference be- 
tween the "p.s.f." and model magnitudes. In addition, we 
trained our photometric-redshift neural network to improve 
the efficiency of star-galaxy separation. We introduced an 
additional parameter (5 sg to the network such that <5 sg = 1 
if the training 2SLAQ object is a galaxy and <5 sg = if the 
object has been spectroscopically identified as a star. The 
network then output values of 8 sg for each photometric ob- 
ject. These values form a continuous distribution between 
and 1 and may be interpreted as a star-galaxy classification 
probability. Catalogue galaxies with low values of <5 sg were 
then excluded from the clustering analysis. We chose the cut 
<5 sg > 0.2, which reduced the level of stellar contamination to 
~ 1.5% with the loss of only ~ 0.1% of the genuine galaxies. 
The final catalogue analyzed in this paper contains 644,903 
entries. In Section 4 we consider the effects on the power 
spectrum of varying the star-galaxy separation procedure. 

2.2 Angular selection function 

We measured the clustering of galaxies selected from the 
SDSS 4th Data Release (DR4) in the North Galactic Plane 
(NGP). We analyzed galaxies in the area bounded by right 
ascensions 110° < a < 270° and declinations — 5° < 5 < 
70°. This area excludes the three SDSS stripes in the South 
Galactic Plane (stripes 76, 82 and 86) which contribute a 
small fraction of the surveyed area and which are widely 
spaced from the rest of the survey region. 

We determined the angular selection function within 
the region (110° < a < 270°, -5° < 8 < 70°) using 
the coverage mask provided by the SDSS team (the files 
atStripeDef .par and tsChunk. dr4 .best .par available via 
the web-page http://www.sdss.org/dr4/coverage). The 
publicly-available information includes the extent of the ob- 
servations in each survey stripe and a list of holes due to 
faults or poor-quality data. This information is given in 
terms of SDSS "survey co-ordinates" (pi and v for each 
stripe) which may be converted to right ascension and decli- 
nation using the relations given on the web-page. This data 
was used to produce our initial window function by assign- 
ing "1" to the surveyed areas and "0" to the unsurveyed 
areas (Figure 1). The sky area encompassed by this window 
function is 5,914 deg 2 , mapping out a total cosmic volume 
2.5 fr~ 3 Gpc 3 across the redshift range 0.4 < z < 0.7. 

There are various effects which may cause the true an- 
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Figure 1. The angular coverage map of the NGP region of the SDSS DR4, which was analyzed in this study. The right ascension and 
declination values are plotted in degrees and an equal-area Aitoff projection is used. The total area of this survey region is 5,914 deg 2 . 
The dashed lines indicate constant Galactic latitudes of 45° (the innermost pair), 30° and 15°. There are several small gaps within the 
observed regions due to faults or poor-quality data. 



gular selection function for the survey data to deviate from 
this binary coverage map (Scranton et al. 2005). These ef- 
fects include: varying completeness in the overlap regions 
between survey stripes, dust extinction or seeing variations, 
incompleteness in the vicinity of very bright stars or galax- 
ies, and systematic effects induced by star-galaxy separa- 
tion. If left uncorrected, these effects may potentially cause 
a systematic error in the measured clustering pattern, par- 
ticularly on large scales. We investigate the importance of 
these effects in Section 4 using an extensive series of tests. 
Our conclusion is that, as far as we can tell from these tests, 
the binary coverage map is an acceptable approximation in 
the galaxy magnitude range 17.5 < idev < 19.8 analyzed in 
this study. 



2.3 Photometric redshifts 

We determined photometric redshift estimates for the LRG 
catalogue using the software package "ANNz" (Collister & La- 
hav 2004). In brief, artificial neural networks are applied to 
parameterize the relation between the galaxy redshift and 
the input information (principally the galaxy photometry, 
but other data such as the angular radius containing a given 
fraction of the galaxy flux can be readily included) . The neu- 
ral network is "trained" using a set of galaxies with known 
spectroscopic redshifts (in this subset of the 2SLAQ 

database) by minimizing a "cost function" (essentially, the 
sum of the squared differences between the photometric and 
spectroscopic redshifts). After each training iteration, the 
cost function is also evaluated for a separate "validation" 



set of spectroscopic redshifts (a second subset of the 2SLAQ 
database) in order to avoid over-fitting to the training set 
given the presence of noise. The final photometric-redshift 
performance is objectively assessed by a third subset of the 
spectroscopic sample called the "evaluation" set. Further de- 
tails on the ANNz software may be found in Collister & Lahav 
(2004). 

For the input photometry we used the de-reddened griz- 
band "model magnitudes" of the SDSS galaxies, which pro- 
vide the most robust colour estimates. We excluded the it- 
band photometry owing to concerns over a time- varying red 
leak in this filter that could potentially introduce systematic 
errors in the photometric redshifts as a function of galaxy 
position on the sky. In fact, irrespective of this concern, the 
u-band does not contribute any measurable benefit to the 
photometric redshift accuracy owing to the relatively low 
signal-to-noise ratio of the it-band photometry. 

The resulting photometric redshifts for the conservative 
version of the catalogue analyzed in this paper, assessed us- 
ing the spectroscopic evaluation set, possess a standard red- 
shift error 



(Sz 



Sz) 



1/2 



0.041 



(1) 



where Sz = z p ho 



and do not suffer from any signif- 



icant bias in the mean difference with respect to the spec- 
troscopic redshifts, Sz < 0.001 (considering galaxies in the 
range 0.4 < z apcc < 0.7, 0.4 < z phot < 0.7 and i dcV < 19.8). 
The performance of our code is hence comparable to the re- 
sults obtained for a similar sample of galaxies by Padman- 
abhan et al. (2005). Figure 2 is a scatter plot of photomet- 



© 0000 RAS, MNRAS 000, 000-000 



Cosmological baryonic and matter densities from SDSS LRGs 5 



cn o 



in 
d 



mean (6z) = 
r.m.s. (6z) = 



-0.0005 
0.0409 



'0.4 



0.5 



0.6 



Spectroscopic redshift 



Figure 2. Comparison between the photometric redshift z p hot 
and spectroscopic redshift z spcc of galaxies in our evaluation set 
(for the conservative version of the catalogue analyzed in this 
paper). The mean and standard deviation of the quantity Sz = 



0.45<z ph <0.50 ,14=0.474 CT=0.035 N=2656 

0.50<z ph <0.55 ,u=0.525 <j=0.042 N=2760 

0.55<z ph <0.60 ^=0.572 a=0.044 N=1854 

0.60<z h <0.65 ,11=0.625 CT=0.053 N= 825 




z phot 



arc indicated. 



Spectroscopic redshift 

Figure 3. Gaussian fits to the spectroscopic redshift distribution 
of Luminous Red Galaxies in each of four photometric redshift 
slices analyzed in this study. The best-fitting Gaussian parameters 
\i and a are indicated, where p(z) oc exp {— [(z — fj,) 2 /2a 2 ]}. N is 
the number of galaxies in the evaluation set used to measure the 
distribution in each slice. 



shift slices of width Az = 0.03. The characteristic patterns 
of large-scale structure may be readily discerned. 



ric and spectroscopic redshifts in the range of interest. The 
principal systematic effect is that galaxies with z spcc < 0.45 
are scattered up in photometric redshift such that there are 
very few galaxies with z p hot < 0.45. The robustness of the 
photometric redshifts worsens in the vicinity of z ~ 0.4 be- 
cause the 4000A break enters a region of poor sensitivity 
between the g and the r bands. In addition, galaxies with 
2 sp cc > 0.65 are typically scattered down in redshift. Photo- 
metric redshift performance improves steadily with bright- 
ening magnitudes (and hence decreasing redshifts) , reaching 
(7 z ~ 0.02 for the brightest galaxies. We refer the reader to 
Collister et al. (2006) for more discussion. We note that the 
small level of systematic photo- z bias evident in Figure 2 is 
not a concern because it can be quantified using the spectro- 
scopic redshift data and folded into our survey simulations. 

Figure 3 is a plot of the distribution of spectroscopic 
redshifts in the evaluation set for galaxies binned in four 
photometric-redshift slices of width Az — 0.05 from z = 
0.45 to z = 0.65. We find that these error distributions are 
well fitted by Gaussian distributions; the best-fitting values 
of the mean and standard deviation for the different slices 
are indicated in the Figure. These Gaussian functions are 
taken as our model for the redshift distribution of galaxies 
in each photo- z slice when analyzing the galaxy clustering 
results. We repeated our cosmological analyses using the raw 
binned redshift distributions in each slice in place of the 
Gaussian fits, and found no significant difference in results. 
The number of galaxies with photometric redshifts outside 
the range 0.45 < z < 0.65 is very small. 

The photometric redshift performance is sufficiently 
good that we can observe the clustering visually. Figure 4 
plots the 2000 most luminous galaxies in the stripe 135° < 
q < 235°, 0° < 5 < 10° in each of a series of narrow red- 



3 MEASUREMENT OF THE ANGULAR 
POWER SPECTRUM 

3.1 Method 

A distribution of galaxies can be related to its angular power 
spectrum Ce in two statistical steps. Firstly, the galaxy den- 
sity field projected onto the sky, a(9,4>), is expanded in 
terms of its spherical harmonic coefficients ae im - 

a(9,<P) = J2 ^ ai, m Y t<m {0A) (2) 

1=0 m=—i 

where Ye tTn axe the usual spherical harmonic functions. Sec- 
ondly, galaxy positions are generated in a Poisson process as 
a (possibly biased) realization of this density field. The an- 
gular power spectrum Ci prescribes the expectation values 
of the spherical harmonic coefficients in the first step of this 
model. It is defined over many realizations of the density 
field (indicated by angled brackets) as: 

< \a e , m \ 2 >= Ci . (3) 

The measurement of the angular power spectrum from a 
set of discrete galaxy positions has been discussed by many 
authors. The available techniques can be crudely divided 
into two principal approaches. Firstly, the power spectrum 
can be derived by direct spherical harmonic estimation, in 
which the harmonic coefficients are calculated by summation 
over the galaxy angular positions (e.g. Peebles 1973; Scharf 
et al. f992; Wright et al. 1994; Wandelt, Hivon & Gorski 
2001). Secondly, a maximum likelihood approach can be uti- 
lized (similar to those commonly employed for the analysis of 
CMB temperature and polarization maps). Examples of the 
application of maximum likelihood methods to galaxy clus- 
tering measurements are the papers by Efstathiou & Moody 
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Figure 4. The locations of the 2000 most luminous galaxies in each of a series of photometric redshift slices of width Az = 0.03 in a 
narrow survey stripe. The characteristic patterns of large-scale structure may be observed. The axes are right ascension and declination 
in decimal degrees. 



(2001) ; Huterer, Knox & Nichol (2001) and Tegmark et al. 

(2002) . 

We used the technique of spherical harmonic estimation 
in the present study. We give a brief summary of the method 
here, referring the reader to Blake, Ferreira & Borrill (2004) 
for more details. The spherical harmonic coefficients of a 
density field for a completely observed sky may be estimated 
by summing over the N galaxy positions (8i, (pi): 



A e ,m = ^2 Ye * 



The estimator for the angular power spectrum is then 



\Al, T , 



N/AQ 



(4) 



(5) 



such that < CfZ > 



Ce ■ Afi is the total survey area, which 
is 4-7T for a complete sky. The discreteness of the distribution 
causes the correction term "— N/AQ". For a given multipole 
£ there are 2^ + 1 different estimators of Ce corresponding 
to m = -i, -t + 1, -1, 0, 1, ...,£. 

For an incomplete sky, these equations must be cor- 
rected for the unsurveyed regions, such that an estimate of 
C e is 



^obs _ \Ai t m — (N/AQ,) It, m y 

W,m — j 

Jl.m 



N 
AO 



(6) 



(e.g. Peebles 1973 equation 50) where 



Ye m dQ 



Je,m — 



We. 



dQ, 



(7) 



(8) 



where the integrals are performed over the survey area AQ,, 
and are determined in our analysis by numerical integra- 
tion. For the SDSS area analyzed in this study (Figure 1), 
Af2/47r = 0.143. 

We determined the observed angular power spectrum 
for the ^th multipole, CJ? bs , by averaging equation 6 over m: 



C? 



Et /-to 
m=-e °i 

2£ + l 



(9) 



The partial sky coverage has the effect of convolving the 
harmonic coefficients such that the measured angular power 
spectrum at multipole I depends on a range of multipoles £' 
of the underlying power spectrum: 



Ce, 



(10) 



The "mixing matrix" Re,e' ma y be determined from the an- 
gular power spectrum We of the survey window function: 
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Multipole I 

Figure 5. The mixing matrix Re e> for multipole I = 200 as a 
function of £' , indicating the range of multipoles correlated by the 
survey window function effects (see equations 10 and 11 for the 
definition of Re e')- 



£ £' £" 




(11) 



where the matrix in equation 11 is a Wigner 3-j symbol and 



We = 



2£ + l 



(12) 



(e.g. Hivon et al. 2002; Deligny et al. 2004). In Figure 5 we 
plot normalized values of the mixing matrix Re,e' evaluated 
using equation 11, corresponding to the SDSS DR4 geometry 
displayed in Figure 1, for multipole t = 200 as a function 
of £' . The convolution function 7?200,f is sharply peaked, 
dropping to 10% of its peak value for an offset A£ = 4 from 
the peak. The shape of this function is largely insensitive 
to the choice of I = 200: to a very good approximation, 
changing this value of £ only causes a translation of the 
function along the £'-axis. 

We averaged the measured angular power spectra in 
multipole bands of width Al = 10. It can be seen from 
the foregoing argument that these bands are approximately 
statistically independent (correlated at a level below 10%) 
and we neglected the residual correlations. We determined 
the angular power spectrum up to a maximum multipole 
^max = 500 (in Section 3.4 we limit the range of multipoles 
to be fitted by our theoretical model in accordance with the 
approximate extent of the linear clustering regime). 



3.2 Error determination 

We compared various different determinations of the sta- 
tistical error a(Ce) in our measurement of the angular 
power spectrum. Firstly, the application of Gaussian statis- 
tics leads to a simple analytical approximation (e.g. Dodel- 
son 2003 p.342 eq.11.27): 



/sky (21- 



Ct + 



N/Afl 



(13) 



where N/Afl is the average source density (in units of sr 1 ) 
and / s ky denotes the fraction of sky covered by the survey. 



The structure of this equation can be readily understood in 
crude terms. The contribution to the error for each value 
of m consists of a cosmic variance term (Ce) and a shot 
noise term [(N/AQ,)^ 1 ]. The averaging over m reduces the 
variance by a factor 2£ + 1. The fraction of sky covered by 
the survey enters as a(Ce) oc 1/ \/ / s k y . (If the observed sky 
area is halved and just one of these sub-areas is analyzed, we 
would expect the standard deviation of the measurement to 
increase by a factor \/2, under the approximation that the 
two sub-areas are independent.) 

We used equation 13 in our cosmological parameter fit- 
ting. We note that the value of Ce in the equation could 
correspond to either the measured power or the model power 
at the multipole in question. In our treatment we used the 
model angular power spectrum for each combination of pa- 
rameters to re-derive the statistical error and assess the 
goodness-of-fit of that set of parameters. 

Estimating the measurement error using equation 13 
involves at least two approximations: 

• The effects of the survey window function are encapsu- 
lated in just one quantity / s k y , thus equation 13 is not likely 
to be valid for multipoles £ corresponding to angular scales 
6 ~ 180° /£ which are similar to or larger than those defin- 
ing the structure of the window function (for which edge 
effects, i.e. the distribution of the survey area, are likely to 
be important). 

• Gaussian statistics are assumed, thus equation 13 is not 
likely to be valid in the regime of non-linear clustering. 

However, it may be hoped that equation 13 is a good ap- 
proximation for multipoles £ corresponding to angular scales 
9 smaller than the typical structure of the window function 
(8 < 2.5° in our case), but large enough that the corre- 
sponding spatial scales lie in the (quasi-) linear regime where 
Gaussian statistics should apply. This turns out to be the 
case for our analysis. 

In order to test the validity of equation 13, we estimated 
the error in the angular power spectrum using two additional 
methods: 

• We generated Monte Carlo realizations of galaxy dis- 
tributions consistent with a model angular power spectrum 
similar to that measured for the real data, using the same 
survey window function as plotted in Figure 1. For each re- 
alization we calculated the density field <t(#,</>) (over a fine 
grid of pixels) by evaluating equation 2 for a set of spherical 
harmonic coefficients corresponding to a Gaussian realiza- 
tion of the underlying Ce spectrum. We then produced a 
catalogue of discrete galaxies as a Poisson realization of this 
density field, and measured the angular power spectrum us- 
ing the same techniques utilized for the real survey data. The 
scatter of the measured power spectra across these Monte 
Carlo realizations will include the full effects of the survey 
window function, unlike equation 13, and this investigation 
is therefore designed to assess the approximation of using 
the single quantity / s k y to represent these effects. 

• Using a "boot-strap" approach, we divided the observed 
survey area into many 10 x 10 deg patches (the choice of 
patch size is a compromise between capturing a sufficient 
number of large-scale power spectrum modes within a sin- 
gle patch, and ensuring there are enough patches across the 
sky to permit reasonable statistics) . We measured the angu- 



© 0000 RAS, MNRAS 000, 000-000 



8 Chris Blake et al. 



Monte— Carlo realizations 
Flat-sky sub-fields 




100 



200 



300 



400 



500 



Multipole I 



Figure 6. The ratio of the error in the angular power spectrum 
measurement to the Gaussian prediction of equation 13 deter- 
mined by two different techniques: firstly by Monte Carlo real- 
izations (the solid line), and secondly by an analysis dividing the 
survey area into many independent patches (the dashed line). 
The Gaussian prediction is a good approximation for the scales 
of interest. The "spikiness" of the dashed line is caused by the 
poor signal-to-noise ratio of the measurements in the 10 X 10 deg 
patches, but the average behaviour is in good agreement with the 
Gaussian prediction. 



lar power spectrum of the galaxies within each patch using 
a flat-sky approximation (i.e., by taking a 2D Fast Fourier 
transform of the overdensity field 8 X and evaluating the 2D 
power spectrum P2{k) = \Sk\ 2 - In the small-angle approxi- 
mation, Ci = P^ii).) The standard deviation of the power 
spectra across these patches is not sensitive to the global sur- 
vey window function, but should trace non-Gaussian effects 
present in the data, thereby assessing the approximation of 
Gaussianity in equation 13. 

In Figure 6 we plot the ratio of the angular power spectrum 
error measured by these two techniques to the prediction 
of equation 13. It can be seen that equation 13 is a good 
approximation for the standard deviation for the scales of 
interest, and it is used in our analysis. Section 3.3.2 presents 
an additional test for the Gaussianity of the clustering fluc- 
tuations. 

3.3 Results 

3.3.1 Angular power spectra 

We measured the angular power spectrum of the LRGs (us- 
ing the estimator defined by equations 6 to 9) in four photo- 
metric redshift slices of width Az = 0.05 between z = 0.45 
and z — 0.65. The results, averaged in multipole bands of 
width A£ = 10, are listed in Table 1, together with the Gaus- 
sian errors of equation 13 (divided by V A£ to account for 
the binning). For the purposes of this Table the measured 
power spectra are used to assign the errors with equation 13. 
When fitting cosmological models to these power spectra in 
Section 3.4, the model power spectra are used to determine 
the Gaussian errors. The surface densities N / AO. of galaxies 
in the four slices are (35.7, 29.1, 18.8, 8.7) deg"" 2 . The results 
are also plotted in the left-hand panels of Figure 10. 



3.3.2 Probability distribution of Ae t , n 

The distribution of values of At.m (the estimated spherical 
harmonic coefficients defined by equation 4) is an interest- 
ing probe of the galaxy distribution (Peebles 1973; Hauser 
& Peebles 1973). Consider first an unclustered distribution 
of points with surface density N/AQ over a full sky. The 
central limit theorem ensures that the real and imaginary 
parts of Ae, m — 2j ^7m(*) are drawn independently from 
Gaussian distributions such that the normalization satis- 
fies < \Ai, m \ 2 >= N/AQ. It is then easy to show that 
x = Ae yYn \ 2 has an exponential probability distribution P(x) 



for m^0: 
P[x) dx = 



exp (—x/a) 



dx 



(14) 



where a — N/AQ. For a partial sky, \A(_ :m \ 2 should be re- 
placed by x — \Ai.m — (N/AQ)Ii,m\ 2 / Je, m - For a clustered 
distribution, a assumes the value (N/AQ) + Ce and equation 
14 holds if each Fourier component has a randomly-assigned 
phase. The validity of equation 14 is therefore a good cross- 
check of the applicability of the linear regime, for which the 
random-phase assumption should apply (Chiang & Coles 
2000). 

Figure 7 plots the distribution of the quantity x for har- 
monic coefficients in the range 1 < m < £ in six different 
multipole bands of width A£ = 10, for the first of the photo- 
metric redshift slices defined above (0.45 < z < 0.5). We also 
display the expected probability distribution for random- 
phase statistics, deriving the quantity a = (N/AQ) + Ci 
from the galaxy surface density together with the measured 
strength of the angular power spectrum in the multipole 
band. We find that the observed distribution of harmonic 
coefficients is a very good match to the random-phase pre- 
diction even for the highest multipoles considered {£ — 300) . 



3.4 Cosmological parameter fits 

3.4-1 Independent redshift slices 

The angular power spectrum Ci is a projection of the spa- 
tial power spectrum of fluctuations at different redshifts z, 
P(k,z), where k is a co-moving wavenumber. The equation 
for the projection is: 



C e = — I Po(k) g e (kf dk 

7T 



(15) 



where P(k,z) = Po{k) D(z) 2 , with D(z) the linear growth 
factor at redshift z. We note that this decomposition of 
P(k, z) is strictly only valid in linear theory, and its applica- 
tion at smaller scales is an approximation. We have assumed 
a scale-independent bias factor b at redshift z for the galaxies 
with respect to the underlying matter fluctuations. 

The kernel ge(k) is given by (neglecting redshift-space 
distortions) : 



9i{k) = / je(u) f(u/k)du. 



(16) 



Here, je(x) is the spherical Bessel function and f(x) depends 
on the radial distribution of the sources as 



f[z(x)]=p(z)D(z)^y 1 



(17) 
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Value of |A. | 2 (corrected for partial sky coverage) 



Figure 7. The distribution of values of \Ag m | 2 (corrected for partial sky coverage) for harmonic coefficients 1 < m < t in six multipolc 
bands of width A£ = 10 for the photometric redshift slice 0.45 < z < 0.5. The dashed lines are the predictions of random-phase statistics 
based on the observed strength of the power spectrum due to clustering and shot noise. The slope of this model changes across the six 
panels due to the varying amplitude of the angular power spectrum. 



where x(z) is the co-moving radial co-ordinate at red- 
shift z, and p(z) is the redshift probability distribution of 
the sources, normalized such that J p(z) dz = 1 (see e.g. 
Huterer, Knox & Nichol 2001; Tegmark et al. 2002; Blake et 
al. 2004). A good approximation for these equations which 
is valid for moderately large ^>50 is: 

d = b 2 J P(k = £/x,z)x(zy 2 p(z) 2 (^j 1 dz. (18) 

The angular power spectrum will be modified by 
redshift- space distortions on large scales. These distortions 
significantly affect the amplitude of the projected power 
spectrum on large scales ^<50, owing to the relative nar- 
rowness of each redshift slice. The amplitude of the redshift- 
space distortions is controlled by a parameter (5{z) ~ 
Qm(z) ' 6 /b(z), where the quantities on the right-hand side 
of the equation are evaluated at the centre of each redshift 
slice of our analysis. The effect is to introduce an addi- 
tional term to the kernel of equation 16 such that it becomes 
gi{k) + gf(k) where: 

9$(k) = £j™j' e (u)f'(u/k)du (19) 

(Fisher, Scharf & Lahav 1994; Padmanabhan et al. 2006 
eq.26). 

We can use these expressions to fit cosmological pa- 
rameters to the angular power spectrum data in each red- 
shift slice. The cosmological parameters determine both the 
power spectrum of fluctuations P(k, z) and the co-moving 



distance x(z) in the above equations. For the redshift dis- 
tribution of the sources p(z) we used the Gaussian func- 
tions fitted to the training set data (see Figure 3). We de- 
rived model spatial power spectra using the "CAMB" software 
package (Lewis, Challinor & Lasenby 2000) which is based 
on CMBFAST (Seljak & Zaldarriaga 1996), including cor- 
rections for non-linear growth of structure using the fitting 
formulae of Smith et al. (2003) ("halofit=l" in CAMB). We 
outputted power spectra at z — 0.55 (the mean redshift 
under analysis) and scaled these spectra to other redshifts 
using the linear growth factor D(z) 2 . 

For low multipoles £, we calculated the model angu- 
lar power spectrum for a given set of cosmological parame- 
ters using the exact expression of equation 15, incorporating 
redshift-space distortions. For reasons of speed, we substi- 
tuted the expression of equation 18 at higher values of £ 
when the approximation became acceptable. Before compar- 
ing each model with observations we convolved it with the 
survey window function using equation 10. Figure 8 displays 
the modifications to the angular power spectrum caused by 
redshift-space distortions and by the survey window func- 
tion for an example case. Figure 9 plots the kernel ge(k) 
appearing in equation 15 (including the redshift-space dis- 
tortion contribution) for various multipoles £. This kernel 
controls which physical scales k (in Fourier space) contribute 
power to a clustering measurement at a particular multipole 
I. It can be seen that the projection at any given multi- 
pole involves a relatively small range Ak of scales owing to 
the narrowness of the photometric-redshift slice Ax in co- 
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Table 1. The angular power spectrum measurements in four photometric rcdshift slices of width Az = 0.05 between z = 0.45 and 
z = 0.65. This table lists the values of 10 5 X Ce and the corresponding Gaussian error in multipole bands of width Al = 10. The errors 
cr(Ci) are determined for the purposes of this Table by substituting the measured values of Ce in equation 13 and not from the best-fitting 
power spectrum model. 
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ordinate space (the approximate range is Ak/k = Ax/x, 
where k = l/x). 

The galaxy clustering pattern is well-known to be sen- 
sitive principally to the quantity f2 m fzo (where Q, m is the 
matter density and Ho is Hubble's constant), and secondar- 
ily to the baryon fraction S7b/f2 m - For our initial investiga- 



tion we fixed the values of the Hubble constant Ho = 75 
km s _1 Mpc -1 and the tilt of the primordial power spec- 
trum n sca iar = 1 (e.g. following Cole et al. 2005; we consider 
variations in Ho and n sca i ar in Section 3.4.2). We then var- 
ied over a grid the matter density f2 m , the baryon fraction 
/b = Ob/fim, the present-day normalization of the power 
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Figure 8. The model angular power spectrum for the photo- 
metric redshift slice 0.45 < z < 0.5 under various assumptions: 
(1) the small-angle approximation (equation 18); (2) the exact 
expression in real space (equation 15 with the kernel of equa- 
tion 16); (3) after the inclusion of redshift-space distortions (with 
the addition of the kernel of equation 19); (4) after convolution 
with the survey window function using equation 10. The assumed 
cosmological parameters are H m = 0.26, /b = 0.16, h = 0.75, 
as = 1, j3 = 0.3. The small-angle approximation is acceptable 
for the multipole range I > 50; at low multipoles the results differ 
very significantly. 
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Figure 9. The kernel gi(k) + g^(k) appearing in equation 15 
for various multipoles I. This plot indicates the range of Fourier 
scales k contributing to angular power at a given multipole I 
(each solid line is for a different value of i increasing from left to 
right). The dashed line indicates the kernel implied by the small- 
angle-approximation of equation 18. The assumed cosmological 
parameters are O m = 0.26, /b = 0.16, h = 0.75, erg = 1, /3 = 0.3. 
We again see that the small-angle approximation is acceptable for 
the multipole range i > 50. 



spectrum as, and the constant galaxy bias factor b. In the 
linear regime as and b axe degenerate, but to allow for any 
contribution from quasi-linear modes we marginalized over 
both separately, using a flat prior 0.7 < as < 1.1. We deter- 
mined the co-moving distance x(z) assuming a flat universe 
with the remainder of the energy density provided by a cos- 
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mological constant. We calculated the value of the redshift- 
distortion parameter at each model grid point by assum- 
ing /3 = fi m (zo)°' 6 /&, where zo is the central redshift of the 
photo- 2i slice. 

We determined the relative likelihood of each model us- 
ing: 

Likelihood oc |M|" 1/2 exp [-{d T ' M' 1 d) /2] (20) 

where d = Ce l0 b s — C^ mo d is the vector of differences between 
the model predictions and the data points in each band of 
multipoles; and M is the covariance matrix, evaluated for 
the particular model at the grid point using equation 13 
(the chi-squared statistic is \ 2 = d T M~ 1 d). Since we are 
neglecting correlations between multipole bands, the covari- 
ance matrix here is diagonal. 

We assumed that power spectrum modes with spatial 
scales k<0.2h Mpc -1 could be modelled via equation 18 
using a scale-independent bias factor b (this is compara- 
ble to the smallest scales fitted in lower-redshift surveys - 
fcmax ~ 0.15/i Mpc - in Percival et al. 2001 and fc max ss 0.2 h 
Mpc -1 in Tegmark et al. 2004a). We estimated the equiv- 
alent maximum multipole of the angular power spectrum 
using the fact that at redshift z, a multipole £ probes spa- 
tial power on scales of roughly k = l/r(z) (see equation 
18 and Figure 9). This corresponds to f max m 300 at our 
mean redshift z ~ 0.55 and hence for each redshift slice we 
fitted our models to only the first 30 multipole bands of 
width A£ = 10. Our placement of the non-linear transition 
at these angular scales is justified in Figure 10 by the small 
differences for I < ^ max between the best-fitting non-linear 
power spectrum in each slice (generated by "halof it = 1" 
in CAMB) and the corresponding linear power spectrum with 
the same cosmological parameters (generated by "halof it 
= 0" in CAMB). 

The best-fitting values of fi m and ft, at each redshift, 
marginalizing over the other three parameters including as 
and b, are listed in Table 2. The results - f2 m w 0.27 ± 0.05 
and /b ~ 0.17 ± 0.06 - are consistent across the redshift 
slices. Our measurement of as is weak because we are re- 
stricting ourselves mainly to the linear clustering regime, 
where as and b are largely degenerate. In Table 2 we quote 
the best-fitting values of b in each redshift slice assuming 
as = 1. We find 1.5 < b < 1.9, with this high bias reflecting 
the fact that LRGs prefer denser environments (e.g. Zehavi 
et al. 2005). We note that (in the linear regime) these bias 
measurements should be interpreted as constraints on the 
product as x b. 

The bias systematically increases with redshift for two 
reasons: 

(i) Galaxies in more distant redshift slices are preferen- 
tially more luminous (owing to the fixed apparent magnitude 
threshold) and hence more strongly clustered (e.g. Norberg 
et al. 2002). In Table 2 we list the average absolute i-band 
magnitude Mi of galaxies in each redshift slice, calculated 
using a K-correction obtained from a standard Luminous 
Red Galaxy template. 

(ii) In standard models of the evolution of galaxy clus- 
tering, the bias factor of a class of galaxies increases with 
redshift in opposition to the decreasing linear growth factor, 
in order to reproduce the observed approximate constancy 
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Figure 10. Each row of this Figure corresponds to a different photometric redshift slice, indicated by the caption. The left-hand column 
plots the angular power spectrum data (binned in multipole bands of width A£ = 10) and the best-fitting model in each case (solid line) 
including both redshift-space distortions and convolution with the survey window function. The dashed line indicates a model with the 
same cosmological parameters but excluding the non-linear correction. The vertical dotted line indicates the range of multipoles we used 
in our fitting, I < ^ max = 300. This corresponds to an estimate of the extent of the linear regime: spatial scales fc<0.2/i Mpc -1 . The 
right-hand column displays 1-cr and 2-cr probability contours in the plane of Sl^/fim and f2 m , marginalizing over as and b. The solid 
circle indicates the best-fitting parameter combination in each case. The cross displays the best-fitting result when the four redshift slices 
are combined with appropriate covariances. 
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Table 2. Best-fitting values of the cosmological parameters for the modelling of the angular power spectra for multipoles I < £ max = 300 
in the four photometric redshift slices. The best- fitting values and errors for f2 m and fb = Qb/Q m are obtained by marginalizing over the 
other three parameters, including erg and b. The values of the Hubble parameter h and the tilt of the primordial power spectrum n sca j ar 
are both held fixed in each analysis at the values given in the Table. The best-fitting values of the linear bias b are quoted assuming 
erg = 1. We also quote the value of the x 2 statistic for the best-fitting model and the number of degrees of freedom (d.o.f.) in each 
case, as well as the average galaxy absolute magnitude M{. The combined measurements of all four slices are obtained by calculating a 
covariance matrix between redshift slices as described in the text. 



Slice 


h (fixed) 


"scalar (fixed) 


Q m 


fb 


b 


Mi - 5 logic, h 


x 2 


d.o.f. 


0.45 < z < 0.5 


0.75 


1.00 


0.28 ±0.05 


0.16 ±0.06 


1.51 ± 0.03 


-22.45 


18.8 


26 


0.5 < z < 0.55 


0.75 


1.00 


0.26 ±0.04 


0.14 ±0.06 


1.68 ± 0.05 


-22.67 


22.4 


26 


0.55 < z < 0.6 


0.75 


1.00 


0.24 ±0.04 


0.18 ±0.06 


1.74 ±0.06 


-22.92 


18.7 


26 


0.6 < z < 0.65 


0.75 


1.00 


0.28 ±0.06 


0.20 ±0.07 


1.95 ± 0.12 


-23.17 


14.8 


26 


Combined 


0.75 


1.00 


0.26 ± 0.031 


0.16 ±0.036 






82.1 


113 




0.7 


1.00 


0.28 ± 0.032 


0.16 ± 0.038 






82.3 


113 




0.75 


0.95 


0.28 ± 0.031 


0.16 ± 0.039 






81.9 


113 



of the small-scale clustering length (e.g. Magliocchetti et al. 
2000; Lahav et al. 2002). 

In Figure 10 we plot the measured angular power spec- 
tra in the four redshift slices together with the cosmological 
parameter fits. Each row corresponds to a different redshift 
slice as indicated. The left-hand column displays the data 
points and the best-fitting models (solid line). The dashed 
lines illustrate the corresponding power spectrum models 
with the same cosmological parameters but omitting the 
non-linear corrections. The difference between these curves 
indicates the importance of the non-linear correction. The 
vertical dotted line indicates the position of the maximum 
multipole we used in our fitting, £ lnB ^ = 300. The right-hand 
column plots 1-a and 2-a probability contours in the plane 
of fb and Qm, marginalizing over erg and 6. 

In all redshift slices the best models produce good fits 
to the shape of the clustering pattern over scales £ < f ma x = 
300, as indicated by the acceptable values of x 2 listed in 
Table 2. In fact, the model generally remains a good match 
to the data for smaller scales £ > 300 which are not used 
in the fitting. At low multipoles there is a visual suggestion 
of a small excess (1-cr) of measured power with respect to 
the model prediction, but this statement lacks statistical 
significance. 



3.4-2 Combined redshift slices 

The cosmological results from different redshift slices are 
not independent, because galaxies sampling the same cos- 
mic variance are scattered between slices by the photomet- 
ric redshift errors. Therefore we cannot combine the results 
from the different redshift slices by simply multiplying to- 
gether the four probability maps in Figure 10. Instead, we 
derived a full covariance matrix between the measurements 
of the angular power spectrum in multipole bands in all four 
redshift slices. 

Since different multipole bands are uncorrelated within 
each redshift slice (to a good approximation) , the covariance 
matrix has a simple structure of non-zero diagonals with 
entries for multipoles £i = £j for redshift slices i and j. 
These covariances are given by 




Matter density Q m 

Figure 11. The probability distribution for the cosmological pa- 
rameters O m and Qb/f2 m combining the four redshift slices us- 
ing a full covariance matrix, and marginalizing over erg and the 
four linear bias parameters. 1-cr and 2-er contours are displayed. 
The solid and dashed contours assume Ho = 75 and 70 km s _1 
Mpc — 1 , respectively. 

< a c\ >-<c\>< c\ >= - 2 cr (21) 

f s ky(2£ + 1) 

using the same Gaussian approximation as equation 13. 
Here, C\ ] is the cross angular power spectrum between the 
two redshift slices. This can be determined for a given cos- 
mological model by evaluating the projection equation 15 
including a kernel gi(k) for each redshift slice: 

C?f = ?hh J P (k)gmgi(k)dk (22) 

where fa; and bj are the linear bias factors for the slices. 

Using this form of covariance matrix, we re-fitted the 
cosmological models as before. We marginalized over a sep- 
arate linear bias parameter for each redshift slice (bi — » 64), 
thus we are fitting seven parameters (f2 m , fb, erg, 61 — ► 64). 
The probability contours for the matter and baryon den- 
sities are displayed in Figure 11. The resulting 1-cr mea- 
surements of the parameters are fi m = 0.26 ± 0.031 and 
fb = 0.16 ± 0.036 (for fixed values of H = 75 km s' 1 
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Figure 12. The cross angular power spectra measured between all combinations of the four photometric redshift slices. The plotted 
models are not fitted to the data points, but are calculated using equation 22 assuming the best-fitting cosmological parameters from 
the analysis of the auto power spectra. 



Mpc -1 and n sca i ar = 1, marginalizing over erg and the four 
bias parameters). This constitutes a significant detection of 
the influence of baryons on the clustering pattern of galax- 
ies. The dashed contours in Figure 11 illustrate the effect of 
reducing the value of Ho to 70 km s _1 Mpc -1 . The resulting 
shift is exactly as expected given that the clustering pattern 
is sensitive to the combination Q m Ho- We therefore express 
our result as tt m h = 0.195 ± 0.023 [where h = H /(100 
km s" 1 Mpc -1 )]. The dotted contours in Figure 11 are ob- 
tained if Ho = 75 km s" 1 Mpc" 1 but the value of n sca i a r 
is reduced to 0.95. In this case we obtain the marginalized 
error fl m h = 0.21 ± 0.023. 



These parameter measurements agree with the latest 
analysis of the Cosmic Microwave Background radiation 
(the 3-year results of the Wilkinson Microwave Anisotropy 
Probe, see Spergel et al. 2006). The measurements from the 
CMB data alone are $l m /i « 0.17 ± 0.03 and flb/n m ~ 
0.18 ± 0.01. Our results are also in concord with those de- 
rived from the latest spectroscopic redshift surveys. For ex- 
ample, Cole et al. (2005) report Q m h = 0.17 ± 0.02 and 
fib/^m = 0.19 ± 0.05 from an analysis of the final 2dF 
Galaxy Redshift Survey, fixing the Hubble parameter h and 
scalar index of primordial fluctuations n sca i ar in a manner 
similar to our analysis. There is tension at the 1-a level be- 
tween the different best- fitting values for S7 m , but we ascribe 
no statistical significance to this. Our parameter constraints 
are also consistent with the combined analysis of the CMB 
data and power spectrum of the spectroscopic sample of 
SDSS LRGs by Huetsi (2006b). 



3.4-3 Cross power spectra 

Equation 22 can be used to evaluate the cross angular power 
spectrum between different redshift slices i and j for a given 
cosmological model. The cross power spectrum can also be 
measured from the data; this is a useful cross-check on the 
consistency of our analysis (i.e. the best-fitting cosmological 
parameters and the redshift distributions within each photo- 
z slice). The estimator is: 



cr = 



2£+l ^ 



(A 



A-i >m 



(23) 



where A\ m and A\ are the spherical harmonic coefficients 
estimated in redshift slices i and j, corrected for partial sky 
effects as described in Section 3.1. By analogy with equation 
13, the variance in this estimator is: 



/sky {21 + 



(24) 



where C\ and C\ are the auto angular power spectra in the 
two redshift slices and Ni and Nj are the numbers of galaxies 
analyzed. 

Figure 12 plots the measured cross power spectra for 
the six possible combinations of four redshift slices, together 
with the model predictions for the best-fitting cosmological 
parameters from the combined analysis of the auto power 
spectra of all the redshift slices. Although the cross power 
spectra have not been used in the cosmological fitting, the 
agreement is generally excellent. The only combination of 
redshift slices which performs poorly is 1 and 4, the most 
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Figure 13. Combination of the angular power spectrum mea- 
surements. The data is first divided by a smooth polynomial fit 
and then scaled along the x-axis in order to match equivalent 
spatial scales projected from different redshifts. The aim is to 
detect "model-independent" evidence for baryon oscillations in 
the clustering pattern. The solid sinusoidal line is a model of the 
result. 



4 TESTS FOR SYSTEMATIC PHOTOMETRIC 
ERRORS 

In this Section we perform an extensive series of tests for 
potential systematic photometric errors that may affect our 
clustering results. Unfortunately, any small photometric er- 
rors will be amplified for a sample of Luminous Red Galaxies 
because these objects are selected from a very steep part of 
the galaxy luminosity function, where a relatively small shift 
in the magnitude completeness threshold (e.g., Am = 0.05) 
will produce a much larger (10s of per cent) change in the 
galaxy surface density. In order to search for such effects, 
we compared the angular power spectrum measured for the 
"default" sample (defined in Section 2) with that obtained 
by restricting or extending the galaxy selection in the fol- 
lowing ways: 

• Exclusion of areas of high dust extinction. 

• Exclusion of areas of poor astronomical seeing. 

• Exclusion of areas in the vicinity of very bright objects. 

• Exclusion of areas lying in the overlap regions between 
survey stripes. 

• Exclusion of areas of low Galactic latitude. 

• Separate analysis of the two largest disconnected survey 
regions in Figure 1. 

• Variations in the star-galaxy separation criteria. 



separated pair. In this case, the cross-correlation amplitude 
is generated by the furthest wings of the redshift probability 
distribution in each slice. In this regime the Gaussian model 
falls off too quickly, hence the model prediction lies below 
the cross power spectrum data points. 



3.5 Searching for baryon oscillations 

We used the angular power spectrum measurements to look 
for "model-independent" evidence of baryon oscillations in 
the clustering pattern. We divided out the overall shape 
of the power spectrum in each case via a smooth polyno- 
mial fit, and searched for significant enhancements or re- 
ductions in power at the expected scales. The baryon oscil- 
lations encode a preferred spatial scale (fixed in co-moving 
co-ordinates) as an approximately sinusoidal modulation of 
power in Fourier space. Therefore, given data with a suffi- 
ciently high signal-to-noise ratio, we should observe the cor- 
responding preferred angular scales moving to larger values 
of £ with increasing redshift, which would constitute an ex- 
cellent "model-independent" test for the existence of baryon 
oscillations. Unfortunately, the current data does not pos- 
sess a sufficiently high signal-to-noise ratio to permit this 
test. We therefore combined the angular power spectrum 
measurements in the different redshift slices, scaling a mea- 
surement at redshift z along the :r-axis by a factor propor- 
tional to l/x(z) in order to match equivalent spatial scales 
projected from different redshifts. One binning of the data 
is plotted in Figure 13, together with a model prediction. 
We find visual suggestions of baryon oscillations, but the 
statistical significance is less than 3-cr. We have not taken 
into account covariances between redshift slices when con- 
structing Figure 13. 



4.1 Dust extinction 

Figure 14 displays the distribution of dust extinction across 
the SDSS NGP (obtained from the public database). Varia- 
tions in dust extinction may imprint systematic changes in 
galaxy density on the sky near the magnitude limit of the 
survey, or if there are systematic errors in the extinction 
map. In this latter case, dust extinction may affect photo- 
metric redshift estimates by reddening galaxy colours. We 
repeated the angular power spectrum measurement exclud- 
ing regions with high dust extinction (> 0.1 mag; 14% of 
the survey area). The results in each of the four redshift 
slices are shown by the solid circles in Figure 17. These data 
points illustrate the number of standard deviations by which 
the angular power spectrum changes from its default value 
given the variation in the selection criteria. Restricting the 
analyzed area to regions of low dust extinction produces no 
significant change in the power spectrum measurement. 

4.2 Astronomical seeing 

Figure 15 displays the distribution of seeing across the SDSS 
NGP (obtained from the public database). Poor seeing af- 
fects the reliability of the galaxy photometry and the effi- 
cacy of the star-galaxy separation. We repeated the angular 
power spectrum measurement excluding regions with poor 
seeing (> 1.5 arcsec; 15% of the survey area). The result is 
shown by the open circles in Figure 17. Again, any change in 
the value of the power spectrum is less than the statistical 
error in the measurement. 

4.3 Incompleteness near bright objects 

In the immediate vicinity of bright stars or galaxies, the 
depth (completeness) of the SDSS imaging survey decreases 
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Figure 14. Variation in dust extinction across the DR4 NGP region. 
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Figure 15. Variation in astronomical seeing across the DR4 NGP region. 
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Figure 16. Location of overlap regions between SDSS stripes in the DR4 NGP region. 
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Figure 17. The dependence of the angular power spectrum measurement in 4 redshift slices on varying dust extinction and astronomical 
seeing. Results from the default catalogue are compared to an analysis restricting the regions analyzed to (1) a maximum dust extinction 
of 0.1 mag (see Figure 14), or (2) a maximum seeing of 1.5 arcsec (see Figure 15). The x-axis is multipole I and the j/-axis displays the 
number of standard deviations of the new measurement from the default result (with ±ltr marked as horizontal dashed lines). 



Variations in angular selection function (2) 
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Figure 18. The dependence of the angular power spectrum measurement in 4 redshift slices on the presence of bright objects and stripe 
overlap regions. Results from the default catalogue are compared to (1) an analysis placing circular masks of radius 1 arcmin around 
all objects with i-band magnitudes brighter than 12, and (2) an analysis excluding overlap regions between stripes (see Figure 16). The 
results are displayed in the same manner as Figure 17. 
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Variations in angular selection function (3) 
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Figure 19. The dependence of the angular power spectrum measurement in 4 redshift slices on Galactic latitude and disconnected region. 
Results from the default catalogue are compared to (1) an analysis excluding regions below Galactic latitude 45°, and (2,3) a separate 
analysis of each of the two largest disconnected regions in the survey window function (see Figure 1). The results are displayed in the 
same manner as Figure 17. 
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Figure 20. The dependence of the angular power spectrum measurement in 4 redshift slices on the star-galaxy separation criteria. 
Results from the default catalogue are compared to analyses changing the value of the co-efficient C in equation 25. The results are 
displayed in the same manner as Figure 17. 
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owing to the difficulty in distinguishing close neighbours to 
these bright objects against the enhanced background. The 
result is a spurious reduction in the number of very close 
galaxy pairs (with separations 6 < 0.02°) at faint magni- 
tude thresholds. This effect is negligible for our analysis ow- 
ing to the large angular scales we are considering (multipoles 
I < 500 of the angular power spectrum) and to our exclusion 
of bright (idcv < 17.5) and faint (idev > 19.8) objects. In 
order to verify this fact, we repeated the angular power spec- 
trum measurement amending the angular selection function 
by placing a series of circular masks of radius 1' around very 
bright SDSS imaging objects with magnitudes i < 12 (re- 
moving 1.4% of the survey area). The result is shown by 
the solid circles in Figure 18, which differ from the default 
measurement by much less than the statistical error. 

4.4 Stripe overlap regions 

The SDSS is performed in scans which result in long data 
stripes of width 2.5° spanning the sky. These stripes over- 
lap towards the edges of the surveyed region, as displayed in 
Figure 16, causing an increase in the effective depth of the 
survey in these areas of overlap (or equivalently, an increase 
in the completeness at a fixed magnitude threshold near the 
survey limit). We repeated the angular power spectrum mea- 
surement excluding these overlap regions (18% of the survey 
area) . The resulting deviation is shown by the open circles in 
Figure 18. Any differences are small, and we conclude that 
our faint magnitude threshold of 19.8 is sufficiently bright 
that differential completeness effects are unimportant in our 
analysis. 

4.5 Galactic latitude 

The residual stellar contamination in our photometric 
database may be a function of Galactic latitude, potentially 
imprinting large-scale systematic density fluctuations. Lines 
of constant Galactic latitude are overlaid on the survey win- 
dow function in Figure 1. We repeated the angular power 
spectrum measurement excluding areas with Galactic lati- 
tudes below 45° (30% of the survey area). The resulting de- 
viation is shown by the solid circles in Figure 19. Again, the 
differences can be largely accounted for by random statisti- 
cal error, and we conclude that any sky variation of stellar 
contamination has a negligible effect on the clustering mea- 
surements. 

4.6 Analysis of disconnected regions 

The survey window function (Figure 1) is mainly comprised 
of two large disconnected regions. These northern and south- 
ern regions, which we will call A and B, comprise 32% and 
67% of the total survey area, respectively. The lack of over- 
lap between these regions implies that relative photometric 
calibration is challenging. It is plausible that some system- 
atic calibration difference between these two regions could 
imprint apparent large-scale fluctuations in power. We re- 
peated the angular power spectrum measurement for each 
region separately. The resulting deviations are shown by the 
open circles and triangles in Figure 19, which reveal no sys- 
tematic difference in results from the joint analysis of the 



combination of regions. We conclude that the accuracy of 
the relative photometric calibration between these two re- 
gions is acceptable. 



4.7 Star-galaxy separation 

The unavoidable presence of a small fraction / of stars in 
our analyzed catalogue affects the clustering measurement. 
The addition of a distribution of uncorrelated stars to our 
clustered galaxies reduces the overall amplitude of the power 
spectrum. In equation 18, the effect is that the redshift prob- 
ability distribution p(z) is now normalized to (1 — /). To 
first order, the clustering amplitude is reduced by a factor 
(1 — /) 2 , which may be conveniently absorbed into the con- 
stant galaxy bias factor fitted to the data in equation 18. 

There are also second order effects. Firstly, problems 
would arise at low Galactic latitudes owing to the varying 
stellar density with position. This is not important for the 
surveyed SDSS regions, as demonstrated above. Secondly, 
the optical quality of the SDSS camera varies a little with 
camera column (i.e., perpendicular to the stripe direction) 
across its 2.5° diameter. This affects our ability to separate 
stars and galaxies as a function of SDSS ^-coordinate (i.e., 
across stripes). Near the edges of the camera, the point- 
spread function is a little broader, and thus faint compact 
galaxies will be preferentially mistaken for stars, inducing a 
fluctuating galaxy surface density with minima at the edges 
of stripes. The angular scale of the resulting fluctuations 
will be similar to the width of the stripes, 2.5°. Unfortu- 
nately, this is comparable to the preferred angular scale of 
the baryon oscillations (~ 5° at z — 0.5). 

In Figure 21 we illustrate this effect by plotting the 
average overdensity of the LRG catalogue as a function of 
pixel position across the field-of-view of the SDSS camera. 
We bin galaxies as a function of the quantity x = 3000 x 
(n ca m — 1) + n co i, where n cam is the SDSS camera column 
containing the object (an integer between 1 and 6), n co i 
is the pixel position of the galaxy centroid in that column 
(a decimal up to 2048), and 3000 is an arbitrary offset to 
separate the camera columns clearly in the Figure. The solid 
points represent the average galaxy overdensity using the 
star-galaxy separation criteria described by Collister et al. 
(2006), 

Ipsf — imodol > C X (21 — Idov) 

0i, de v > 0.2" (25) 

where the coefficient C = 0.2. In this analysis we do not 
include any additional neural network star-galaxy separa- 
tion using ANNz (see Section 2). For the reasons discussed 
above, the resulting LRG overdensity peaks in column 3 and 
has minima at the camera edges (columns 1 and 6). In or- 
der to correct for this systematic imprint, we tried varying 
the coefficient C in equation 25 to trace the quality of the 
star-galaxy separation, such that C assumed different values 
(0.1, 0.2, 0.25, 0.2, 0.2, 0.16) in the six camera columns. These 
modified values of C are chosen to ensure that the average 
galaxy overdensity is constant across the camera field (the 
open circles in Figure 21). 

We then repeated the angular power spectrum measure- 
ment comparing our catalogue with default star-galaxy sep- 
aration (C = 0.2) to the model where C varies with cam- 
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Figure 21. The average overdensity of LRGs as a function of 
pixel position across the SDSS camera (perpendicular to the di- 
rection of the camera columns). This position is quantified by the 
value of3000x(n cam — l)"t"^col j where 7i C am is the camera column, 
n co i is the pixel position of the galaxy centroid in that column, 
and 3000 is an arbitrary offset to separate the camera columns 
clearly in the Figure. Owing to the varying optical quality of the 
camera across the field-of-view, if constant star-galaxy separation 
criteria are applied then the resulting overdensity profile will ex- 
hibit a peak at the centre of the field with minima at the edges 
(as shown by the solid circles). In order to correct for this effect, 
we varied the aggressiveness of the star-galaxy separation as a 
function of camera column as described in the text, removing the 
systematic density variation (as shown by the open circles). 



era column. We also considered alternative constant values 
C = 0.1 and C — 0.3. The results are displayed in Figure 
20. Our first conclusion is that the adoption of the cam- 
era column-dependent star-galaxy separation is of negligible 
benefit. The imprint of the varying overdensity in Figure 21 
in fact has no measurable effect on the power spectrum. We 
therefore retain our default criteria of a constant coefficient 
C — 0.2. The effect of changing this value is in reasonable 
accord with our model of a constant factor in power, which 
may be safely absorbed into the fitted galaxy bias factor. 
The possible exception to this conclusion is the lowest mul- 
tipole bands in the furthest redshift slice 0.6 < z < 0.65 
(which is expected to contain the most significant stellar con- 
tamination). In this case, more liberal star-galaxy separation 
actually increases the observed power on the largest scales. 
In fact, the lowest multipole band for this redshift slice is 
observed to possess marginally significant excess power in 
Figure 10, thus this one data point may be influenced by 
stellar contamination. 



5 MEASUREMENT OF THE SPATIAL 
POWER SPECTRUM 

In this Section we determine the spatial galaxy power spec- 
trum by a direct three-dimensional Fourier transform of the 
gridded photo-z catalogue. We believe that this is a novel 
technique which has not been previously applied to galaxy 
data, although Seo & Eisenstein (2003) and Blake & Bridle 
(2005) give theoretical discussions. 



5.1 Method 

The photo- z error distribution produces a convolution of the 
underlying density field in the radial direction, and hence 
(in the first approximation) a multiplicative radial damping 
of the measured three-dimensional power in Fourier space. 
In principle this damping factor can be calculated using the 
known photo- z error properties and then divided out, restor- 
ing the underlying unconvolved power spectrum. We must 
implement a flat-sky approximation such that the photo-z 
smearing only occurs along one axis. 

For example, in the ideal case that the radial co-moving 
co-ordinate x of each galaxy is smeared by an amount 8x 
sampled from a Gaussian distribution 



f(Sx) oc exp 



1 fSxy 

2 \V X ) 



(26) 



then the resulting power spectrum signal will be damped in 
the radial direction in accordance with: 



P(k x 



) P(k x 



ky , k- 



) x exp [-(k x UxY 



(27) 



Thus Fourier modes with high values of \k x \ 2> l/cr x con- 
tribute only noise to the measurement. However, if we only 
retain power-spectrum modes with \k x \ < k x ,max ~ 1/fx, 
then we can deconvolve the measured 3D power spectrum 
by dividing by the function exp [—(k x a x ) 2 ] prior to combin- 
ing the Fourier modes in bins of fixed k — y/k x + ky + k% 
(where k x < t im „). In practice the photo-z error distri- 
bution is not a constant Gaussian function, and we calcu- 
lated the deconvolution function for our sample using Monte 
Carlo simulations rather than an analytic approximation, as 
described below. 

We initially analyzed the galaxies in photometric red- 
shift slices of width Az — 0.05, as for the angular power 
spectrum measurement. We implemented the flat-sky ap- 
proximation by dividing the survey area into patches of 
width 20 x 20 deg. These patches have sufficient area to 
contain a reasonable number of large-scale power spectrum 
modes with k < 0.2 h Mpc -1 , and yet are small enough that 
the photo-z smearing may be considered to act along one 
axis of a survey cuboid (which we take as the i-axis). We 
only analyzed patches for which the survey angular selec- 
tion function covered over 50% of the area of the patch. The 
equations mapping the galaxy right ascensions and declina- 
tions (a, S) onto spatial positions (X, Y, Z) in the cuboid 
are 



X 

Y 

Z 



dx 
dz 



(z ) x (z - z ) 



x(zo) x (a — ao) cos<5o 
x{zo) x (5 - S ) 



(28) 



where (cvo,<5o) are the central angular co-ordinates of each 
patch, z is the galaxy redshift, zo is the central redshift of 
each slice, and x(z) is the co- moving radial distance of the 
survey redshift. Thus we must assume values for the cosmo- 
logical parameters to construct this mapping, as discussed 
below. 

We then measured the 3D power spectrum of galaxies 
in each cuboid in Fourier bins of width AA: = 0.01 h Mpc - , 
correcting for survey window function effects and assuming a 
radial selection function which we directly fitted to the sur- 
vey data. We only included a Fourier mode (k x ,k y ,k z ) in 
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the analysis if the radial wavescale \k x \ < /c x ,max = 1.5/a x , 
where a x is the standard deviation of the photo-z errors 
within the slice in co-moving co-ordinates. For redshift slices 
of thickness Az — 0.05 this restricts us to modes with 
k x = 0, because the spacing of the modes in radial fc-space is 
Ak x = 2it/L x > fcx.max, where L x is the width of the slice in 
co-moving co-ordinates (we also calculate the spatial power 
spectrum for a wider redshift slice with Az = 0.15). The 
final power spectrum was determined by averaging the mea- 
surements in the separate patches. The error in the power 
spectrum was taken to be the standard deviation of the mea- 
surements in the patches. 

The calculation of the spatial power spectrum P(k) in 
this manner depends upon the choice of cosmological pa- 
rameters (which specifies the cosmological distances x(z) in 
equation 28). Therefore, the determination of best-fitting 
cosmological parameters using this method must proceed in 
a more complex fashion than the angular power spectrum 
fits in Section 3. In a measurement of P(k), both the data 
and the model depend on the trial cosmological parameters. 
In contrast, the measured angular power spectrum does not 
depend on the cosmological parameters and may be left fixed 
during the model-fitting process. 

However, there is an important advantage in perform- 
ing a direct measurement of P(k) from a photo-z catalogue: 
we avoid unnecessary radial binning of the data. The pro- 
jection of P(k) to the angular power spectrum Ce (equation 
18) corresponds to a convolution of P(k) in Fourier space 
in which the contrast of any features (such as the baryon 
oscillations) will be reduced. But direct manipulation of the 
3D power spectrum involves only a multiplicative damping 
of power (equation 27) which does not diminish the contrast 
of features. Hence the direct power spectrum technique pre- 
sented in this Section may yield the most significant detec- 
tion of baryon oscillations. In this initial study we do not in 
fact fit cosmological parameters to the P(k) measurements, 
but instead assume the best-fitting set of parameters from 
the angular power spectrum analysis (Table 2) in order to 
demonstrate the consistency of the results and to search for 
the existence of baryon oscillations. 

The spatial power spectrum measured by this direct 
Fourier calculation must be carefully normalized when the 
condition a x <C L x is not satisfied, i.e., when the photo-z 
error is comparable to the width of the survey slice (which 
is the case here). For a spectroscopic survey of volume V, 
the dimensionless power spectrum of the Fourier amplitudes, 
<5fc| 2 , must be multiplied by V in order to obtain a measured 
quantity (e.g. with units of h~ 3 Mpc 3 ) which is independent 
of the volume of the survey and may be compared with 
theoretical predictions. Suppose this spectroscopic survey is 
now convolved with a photo-z error function, and a radial 
slice of the resulting photo-z catalogue is analyzed, with vol- 
ume Vsiico- The correct normalization factor for comparison 
of the measured power spectrum with theory is not in fact 
Vsiico, because the galaxy distribution within the slice has 
been smoothed over a wider volume owing to the scattering 
of objects into the slice by the photo-z error distribution. 
In the ideal case of a "top-hat" radial selection function for 
the photo-z slice (of width L x ) and the Gaussian photo-z er- 
ror function of equation 26, the multiplicative normalization 
correction is given by 



2tt 



du exp 



( u V" 




sin (it/2) 


\L x /a x J 




[ («/2) J 



(29) 



such that the effective normalization volume for the power 
spectrum measurement is V B u cc /N. Equation 29 has the lim- 
its TV — > 1 for a x <C L x (i.e., a spectroscopic survey) and 
TV — > L x /(2y^ra x ) for L x < a x . 

Moreover, the smoothing of the density distribution 
from outside the photo-z slice affects the shot noise cor- 
rection term P s hot subtracted from the raw power spec- 
trum measurement, |<5fc| 2 . The shot noise correction is not 
simply given by the volume density n of analyzed galax- 
ies within the slice, P s hot = 1/n. Thus the measured power 
spectrum P mC as must be adjusted by both additive (shot 
noise) and multiplicative (normalization) corrections in or- 
der to achieve a reliable estimate of the underlying power 
spectrum P tr uc: 

Ptruc (fc) = A(k) P meas (fe) + B(k) Pshot (30) 

where A(k) and B(k) are roughly constant with k = 
■\/fc| + ky + k\ , although change rapidly with k x (< 
fez, max <C k) owing to the exponential damping factor in 
equation 27. 

We determined the two correction functions A(k) and 
B(k) in equation 30 using Monte Carlo realizations. We 
first generated clustered distributions of galaxies in a three- 
dimensional survey cone, using a model spatial power spec- 
trum corresponding to the best-fitting cosmological parame- 
ters. These Gaussian power-spectrum realizations were cre- 
ated using the method described in Glazebrook & Blake 
(2005). The radial selection function before the application 
of photometric redshift errors was fixed by the spectroscopic 
redshift distribution of the 2SLAQ training set (our best es- 
timate of the unconvolved redshift distribution). The SDSS 
angular selection function of Figure 1 was then imposed. 
We then applied a photometric redshift convolution func- 
tion, again determined using the training set, such that a 
different photo-z error distribution is implemented at each 
spectroscopic redshift, to ensure maximum consistency with 
the statistics of the real catalogue. The appropriate photo- 
z slices were then isolated and the spatial power spectra 
were measured with the same methods as applied to the 
real data. In order to complete the calculation, a second set 
of Monte Carlo realizations was generated using an identi- 
cal method, except that the original spectroscopic catalogue 
was unclustered (i.e., Ptruc = in equation 30) and thus 
the measured power spectrum purely probes shot noise ef- 
fects. The two measured power spectra, of the clustered and 
unclustered simulations, in conjunction with the known in- 
put power spectrum, allow us to deduce the functions A(k) 
and B(k) in equation 30. These correction factors are then 
applied to our power spectrum measurements from the real 
data. 



5.2 Results 

We display our spatial power spectrum measurements in the 
four redshift slices in Figure 22, together with the models 
corresponding to the best-fitting cosmological parameters 
from the angular power spectrum analysis (Table 2). The 
agreement is good, validating our analysis techniques. There 
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Figure 22. The spatial power spectrum of the photo-z catalogue in each of four redshift slices of width Az = 0.05, measured using direct 
Fourier techniques in a series of 20 X 20 deg patches. There is no data point in the first bin (at k = 0.005 h Mpc -1 ) because this bin 
contains no Fourier modes in the patches. The solid lines are not a fit to this data, but are the model power spectra corresponding to the 
best-fitting cosmological parameters (C m , /b, o"8> b) for the angular power spectrum analysis in each redshift slice. The good agreement 
between the data and model validates our analysis techniques. 



is evidence of excess power in the largest-scale bin, possi- 
bly due to the failure of the small-angle approximation of 
equation 28 on these scales. We note that redshift-space dis- 
tortions do not affect the power spectrum models plotted 
in Figure 22 because we are only analyzing power spectrum 
modes with k x — 0. The power spectrum boost due to large- 
scale coherent infall depends on the angle to the line-of-sight 
9 of a given Fourier mode as (1 + (3 cos 2 9), thus modes with 
k x ~ (i.e. 9 « 90°) contain little imprint of redshift-space 
effects. 

In Figure 23 we plot the results of analyzing a wider red- 
shift slice of the survey, 0.475 < z < 0.625, such that Fourier 
modes with k x ^ are being utilized. Again, agreement with 
the model spatial power spectrum is very good. The shape 
of our recovered power spectra also agree well with the mea- 
surement from the spectroscopic sample of SDSS LRGs at 
lower redshift by Huetsi (2006a). 

Unfortunately the current size of the survey is not suffi- 
cient to yield convincing "model-independent" evidence for 
baryon oscillations. However, a modest expansion of the sur- 
vey volume (e.g. the final SDSS Data Release) would clearly 
create the potential to produce such a detection (as mod- 
elled by Blake & Bridle 2005). 



6 CONCLUSIONS 

We have quantified the patterns of large-scale structure 
in a photometric-redshift catalogue of intermediate-redshift 



0.475 < z < 0.625 




0.1 0.2 



k / h Mpc 1 

Figure 23. The spatial power spectrum of the photo-z catalogue 
in a wide redshift slice Az = 0.15 such that Fourier modes with 
k x ^ are utilized in the analysis. The agreement is good with the 
model power spectrum corresponding to the best-fitting cosmo- 
logical parameters for the angular clustering measurements (the 
solid line). 

(0.4 < z < 0.7) Luminous Red Galaxies from the SDSS 
4th Data Release. This catalogue, which we have named 
MegaZ-LRG, contains > 10 6 photometric redshifts derived 
with ANNz, an Artificial Neural Network code (Collister et 
al. 2006). In this study we analyze a "conservative" version 



© 0000 RAS, MNRAS 000, 000-000 



Cosmological baryonic and matter densities from SDSS LRGs 23 



of the catalogue, which includes 6 x 10 5 objects over 5,914 
deg 2 . The main conclusions of our study are: 

• We measured angular power spectra in a series of photo- 
z slices of width Az — 0.05 from z = 0.45 to z = 0.65. 
The errors in the Ce spectra were tested with Monte Carlo 
realizations and with a "bootstrap" analysis in sub-fields. 
A simple Gaussian approximation was determined to be an 
accurate approximation to the error in Ce over the scales of 
interest. 

• The angular power spectra are well-fitted by model 
spatial power spectra with consistent cosmological param- 
eters n m « 0.27, fb « 0.17, H = 75 km s _1 Mpc -1 and 
^scalar = 1 over the series of redshift slices. 

• Combining the results for the different redshift slices 
with appropriate covariances results in a measurement of 
n m h = 0.195 ± 0.023 and fl b /Q m = 0.16 ± 0.036 (for fixed 
values of Ho and n sca iar , marginalizing over as and b) . These 
measurements agree with the latest analysis of the Cosmic 
Microwave Background radiation (the 3-year results of the 
Wilkinson Microwave Anisotropy Probe, see Spergel et al. 
2006). The measurements from the CMB data alone are 
O m /i = 0.17 ± 0.03 and n b /tt m = 0.18 ± 0.01. 

• These results agree very well with a parallel investi- 
gation of a similar dataset carried out by Padmanabhan 
et al. (2006) using a different photometric-redshift estima- 
tion method and power spectrum analysis procedure. Pad- 
mananhan et al. quote cosmological parameter fits flmh = 
0.21 ± 0.021 and Q b /n m = 0.18 ± 0.04, in both cases within 
1-a of our results. 

• We detect visual suggestions of baryon oscillations with 
a statistical significance of less than 3-cr. 

• We performed a direct re-construction of the spatial 
power spectrum using a Fourier analysis which only retains 
large-scale radial components. We demonstrated that this 
method produces results which are consistent with the best- 
fitting model for the angular power spectrum analysis. 

• We searched for evidence of systematic photometric er- 
rors in the catalogue using a variety of tests, and were unable 
to identify any significant effects. 

• We measure a hint of excess power relative to the best- 
fitting cosmological model on the largest scales (the lowest 
multipole bands in the four redshift slices in Figure 10). If 
confirmed, this excess power has a range of possible causes: 
(1) residual systematic errors; (2) cosmic variance; (3) large- 
scale galaxy biasing mechanisms; (4) new early-Universe 
physics. 

The cosmological parameter measurements we obtain 
from the large-scale clustering of this photometric redshift 
catalogue have a similar precision to those derived from the 
latest spectroscopic redshift surveys, despite the weaker red- 
shift information. For example, Cole et al. (2005) report 
Cl m h = 0.17 ±0.02 and Q b /Cl m = 0.19 ±0.05 from an analy- 
sis of the final 2dF Galaxy Redshift Survey, fixing the Hub- 
ble parameter h and scalar index of primordial fluctuations 
^scalar in a manner similar to our analysis. We therefore con- 
clude that photometric redshift surveys have become com- 
petitive with spectroscopic surveys for the measurement of 
cosmological parameters using the large-scale clustering pat- 
tern in the simple "vanilla" model. 
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